Possible isotope effect on the resonance peak formation in high-T c cuprates 
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Within effective t — J Hamiltonian we analyze the influence of electronic correlations and electron- 
phonon interaction on the dynamical spin susceptibility in layered cuprates. We find an isotope 
effect on the resonance peak in the magnetic spin susceptibility, 1m x(q> w )i seen by inelastic neutron 
scattering, ft results from both the electron-phonon coupling and the electronic correlation effects 
taken into account beyond random phase approximation(RPA) scheme. We find at optimal doping 
the isotope coeffiecient a res m 0.4 which can be further tested experimentally. 
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An understanding of the elementary and the spin exci- 
tations in high-T c cuprates is of central significance. For 
example, it is known that the Cooper-pairing scenario 
via the exchange of antiferromagnetic spin fluctuations 
was quite successful in explaining the various features of 
superconductivity in hole-doped cuprates such as d x i_ y i- 
wave symmetry of the superconducting order parameter 
and its feedback on the elementary and spin excitatitions 
[1] . Most importanly, in this scenario the dynamical spin 
susceptibility, x(q, w), controls mainly the superconduct- 
ing and normal state properties of the layered cuprates 
[1] . One of the key experimental fact in the phenomenol- 
ogy of high-T c cuprates is the occurrence of a so-called 
resonance peak in the inelastic neutron scattcring(INS) 
experiments [2,3]. It occurs below T c in the dynami- 
cal spin susceptibility, %(q, w), at the antiferromagnetic 
wave vector Q = (jr, it) and u> ks uj res which is of the or- 
der of 40meV in the optimally doped cuprates. Its feed- 
back in various electronic properties like optical conduc- 
tivity, Raman response function, and elementary excita- 
tions has been observed experimentally by various tech- 
niques [1]. Furthermore, its successful explanation within 
spin-fluctuation-mediated Cooper-pairing together with 
c? 2 .2_j / 2-wave symmetry of the superconducting order pa- 
rameter favors this scenario as a basic one for supercon- 
ductivity in the cuprates. On the other hand, recent ex- 
periments indicate that also electron-phonon interaction 
influences strongly their behavior [4-8] . In particular, the 
observation of the relatively large isotope effect in vari- 
ous characteristics of cuprates like penetration depth [5] , 
'kink'-structure seen by ARPES [9] still raises a question: 
what is the role of phonons in determining the supercon- 
ducting properties of cuprates? 

Here, we derive an effective t — J Hamiltonian where 
both the hopping integral, t, and the superexchange in- 
teraction between neighboring spins, J, are renormalized 
by phonons. We analyze the influence of the electron- 
phonon interaction on the dynamical spin susceptibility 
in layered cuprates. In particular, we find an isotope 
effect on the resonance peak in the magnetic spin sus- 
ceptibility, Imx(q, ui). It results from both the electron- 



phonon coupling and the electronic correlation effects 
taken into account beyond random phase approxima- 
tion(RPA) scheme. We show that even if the supercon- 
ductivity is driven by the magnetic exchange the char- 
acteristic energy features of cuprates can be significantly 
renormalized by the strong electron-phonon interaction. 

Effective Hamiltonian We start from the atomic limit 
of the three-band p — d Hamiltonian 

H = ^2 e dd\ a d la + ^2 e pPiaP^ + ^2 U d n d ^n dl 

+ U P n Pl n Pl + ^q/q/q U) 



where e d and e p are the on-site energies of the copper and 
the oxygen holes, n da = d\ a d irT and n pa = p\ a Pi a are the 
copper 3d and oxygen 2p hole densities for site i, respec- 
tively. U d and Up refer to the on-site copper and oxygen 
Coulomb repulsion, respectively, /t denotes the phonon 
creation operator and htUq is a phonon energy dispersion. 
The hopping term between copper and oxygen 



Hi 



tpd (dtvPbcr + h.C.) 



and the electron-phonon interaction 
Hx=J2 9^ (/-q + ft 



l — d.p 



(2) 



(3) 



we consider as a perturbation. Here, t pd is a hopping 
term between copper and oxygen, gi is a electron-phonon 
coupling strength at the site /. This notation is similar to 
the simplified Holstein model where the migrating charge 
interacts locally with breathing phonon modes forming 
electron-vibrational states. 

To derive an effective t — J Hamiltonian we employ the 



canonical Schricffcr-Wolf-likc transformations 



,-s 



[10,11]. The matrix of the unitary transformation for the 
initial Hamiltonian is found by excluding the odd terms 
with respect to the hopping integral with an accuracy 
up to the sixth order perturbation theory. Then the S- 
operator consists of the sum of five terms. Each of them 
is determined by the following iteration procedure 



1 



[#oSi] - -H 2 , [H S 2 } = - [ffi5i] , 
[if S 3 ] = - [ffiS 2 ] - ^ [[H2S1] 5i] , 

[ff S 4 ] = - [#iS 3 ] - 5 [[H2S1] S 2 ] - I [[H 2 5 2 ] Si] 

[if S 5 ] = - [Hi5 4 ] - ^ [[H2S1] S 3 ] - i [[H 2 5 3 ] 5i] 

4[[ ff 25 2 ]S 2 ] + ^[[[[ff 2 Si]Si]5i]Si] 



(4) 



The calculations are straightforward and their details will 
be given elsewhere. Note, in the second order perturba- 
tion the effective hopping integral, tij, appears. It is 
further renormalized by the electron-phonon interaction 
in the fourth order term where we introduce the average 
over the phonons. Similarly, the superexchange interac- 
tion occurs in the fourth order perturbation theory and 
its renormalization takes place in the sixth order term. 
Finally, the relevant effective Hamiltonian is given by 



nifij 



(5) 



i>j 



Note, in general case the effective Hamiltonian contains 
also the Coulomb interaction between doped holes and 
the interaction of quasiparticles via the phonon field. We 
dropped these terms here, because they do not contribute 
directly to the spin susceptibility. In Eq.(5) we use the 
projecting Hubbard-like operators ^f' 13 =| i,a >< i,f3 | 
in order to satisfy no double occupancy constraint. The 
index pd corresponds to a Zhang-Rice singlet formation 
with one hole placed on the copper site whereas the sec- 
ond hole is distributed on the neighboring oxygen sites 
[12]. Note, Uj = i9.e-T ls i/fi<''i* where t% is the effective 
hopping integral without taking into account electron- 
phonon interaction, Ei = (g*) 2 /hu>* is the so-called po- 
laron stabilization energy of the copper-oxygen singlet 
state and < 7 < 1. Note, from the experimental data 
[13] the whole exponential factor was estimated to be 
7-E* /huj* w 0.92 around the optimal doping and its value 
is increasing upon decreasing doping. We further assume 
that the lifetime of the Zhang- Rice singlet is much larger 
than the relaxation time of the local deformations. 

Similarly the superexchange interaction between near- 
est copper spins is given by 



J = Jail 



3h 



+E d 0Jd coth 



^ P coth fe) 



2k B T 



(6) 



where A pd = e p — e d + U p — U d is the energy transfer from 
copper to oxygen and is known to be of the order of 1.5eV 
in the cuprates. Here, Jo is the superexchange interaction 
of copper spins via the intermediate oxygen atom in the 
absence of phonons. We took the phonon frequency of 



the order of uj* w uj p w u>d =50mcV which roughly cor- 
responds to the energy of the longitudinal optical(LO) 
bond stretching phonon mode in cuprates. According to 
the recent experiments [6,7] it may play an essential role 
in the physics of cuprates. The so-called polaronic sta- 
bilization energy E* « E p « Ed was estimated of the 
order of 0.5eV in accordance with the measurements of 
the isotope effect in cuprates [14]. 

Dynamical spin susceptibility To derive the dynamical 
spin susceptibility in the superconducting state we use 
the method suggested by Hubbard and Jain [15] that 
allows to take into account strong electronic correlations. 
First we add the external magnetic field applied along c- 
axis into the effective Hamiltonian 



Hi = Re^/ l _ q e^ t - qRi ) + h^ 1 ^'^' 



(7) 



Then we write an equation of motion for the 'J op- 
erators using the Roth-type of the decoupling scheme 
[16] and expanding the P° d = {*J' pd *f < T j = i±*i + 

crRe£ q [s^ q e- 4 (q R j-^) + ^ e <(qR-j -<"*)] U p to the first 
order in S q = x zz (q, w)/i q . In particular, 



Q^pd 



~2 tk -y Sq ~i 



+ 



jTf—VtPd—ito/t 
w k-q c 



<y,pd iu)t 



*k+q-e 



and the similar expression occurs for Here, Ak = 

Aa (cos k x — cos k y ) is e^^-wave superconducting gap, 
Jq = J (cosk x + cos k y ) is the Fourier transform of the 
superexchange interaction on a square lattice. 

The expression for the longitudinal component of the 
dynamical spin susceptibility can be obtained from the 
relation 



(8) 



pd,1 ,t,T ,pd\ 



pd,l,j,l,pd\ 



\ x 1 







(9) 



and using the Bogolyubov-like transformations to the 
new quasiparticle states 



A -k 



Here, 



1 + 



"k* 



pd,a 



k 

a-,pd 



(10) 



and v 2 . = § ( 1 - ) arc the 



"k — 2 ^ 1 

Bogolyubov coefficients, n is a chemical potential, and 
Ek = yj (ik — A*) 2 + A 2 , is the energy dispersion in the 
superconducting state. Substituting Eq.(10) in Eq.(9) 
and using the equations of motion (8) one obtains the 
expression for the dynamical spin susceptibility in the 
form 



x(q,w) 



Xo(q,w) 



JqXo(q, w) + n(q, lu) + Z(q, uj) 



(11) 
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This is a central result of our paper. Here, %o(q, is 
the usual BCS-like Lindhard response function, II(q, u) 
and Z(q,tj) result from the strong electronic correlation 
effects. In the normal state the expression for IT(q, ui) 
has been obtained by Hubbard and Jain [15]. In the 
superconducting state it is given by 



n(q,*) = ^E 



"kWk+q (WkWk+q + ^k^k+q) 



+Vk"k+q (Vk«k+q + WkWk+q) 
+Mk«k+q («k«k+q - Mk+qWk) 
+Mk+qfk (fkWk+q - Mk^k+q) 



^k./k — ^k+q/k+q 

w + i()+ + E k - £ k+q 

*k(l - /k) - ^k+q(l - /k+q) 



+ i0+ - E k + £ k+q 
frc/k - £k+q(l ~ /k+q) 

w + i0+ + E k + £ k+q 

*k(l - /k) - tk+q/k+q 



w + i0+ - E k - £ k+q 
The function Z(q,u) is written as follows 



(12) 



Z(q,w) 



AT ^ 



+ i0+ 



(13) 



Here, /k is the Fermi distribution function, e k ^ = ^^ik, 
e k = Ppdtk is the energy dispersion in the normal 
state, and t k = 2t(cosfc x + cosk y ) + it' cosk x cosk y + 
2t"(cos2k x + cos2k y ) is the Fourier transform of the 
hopping integral on a square lattice including nearest, 
next- and next-next-nearest neighbor hopping, respec- 
tively. The origin of the terms II(q, ui) and Z(q, ui) re- 
lates to the no double occupancy constraint. In partic- 
ular, for the Coulomb repulsion U = oo and J = the 
dynamical spin susceptibility does not reduce to the stan- 
dard Lindhard response function but is renormalized by 
the electronic correlation effects [17]. For the A k = 
Eq.(ll) agrees with the normal state result for the dy- 
namical spin susceptibility [15,18,19]. 

Results and Discussion Inelastic neutron scattering 
(INS) measurements probe directly the imaginary part 
of the dynamical spin susceptibility. Therefore, it is of 
interest to analyze the role played by the electronic cor- 
relations on the 'resonance' peak formation seen by INS 
[3] . Using various approaches this feature was well under- 
stood mainly as a result of the spin density wave (SDW) 
collective mode formation at i.e. when the 

denominator of the spin susceptibility at the antiferro- 
magnetic wave vector Q is close to zero [20,21]. 

In Fig. 1 we show the results of our calculations for the 
Im x(Qj w ) as a function of frequency in the normal and 
the superconducting state. Here, we use t — 200meV, 
t' = — O.lt, and t" = 0.02t at optimal doping. Clearly in 
the normal state the spin fluctuation spectrum is char- 
acterized by a broad feature which starts around lOmcV 
and extends up to a higher frequencies. In the super- 
conducting state it strongly renormalizes due to a pres- 
ence of the d x 2_ y 2-w&ve gap (A k = ^(cosfc^ — cosk y )) 
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FIG. 1. Calculated imaginary part of the dynamical spin 
susceptibility Imx(Q, ui) using Eq.(ll) and J=0.3t in the nor- 
mal and superconducting state at optimal doping. Here, 
we use the superconducting gap Ao = 0.14t (28 meV) and 
T c w 0.04i (90 K) for optimal doping from our earlier cal- 
culations of the mean-field phase diagram [22] . To illustrate 
the role of Z(q, ui) we also show the results for Z(q, ui) = 
(dashed curve). Note, the damping was chosen T = 1.5meV. 



that leads to a resonance peak formation similar to the 
RPA result [20,21]. However, due to a strong frequency 
dependence of Z(Q,w) and n(Q,w) the spectral weight 
of the resonance peak is redistributed away from (it, it) 
(see the dashed curve for comparison) leading to a well- 
pronounced dispersion of the latter. This is illustrated 
in Fig. 2(a) where we show the calculated frequency and 
momentum dependence of Im x(<l> w ) away from the anti- 
ferromagnetic wave vector Q = (tt,tt). The dispersion of 
the resonance excitations is clearly visible and is shown 
in Fig. 2(b) as a function of q x (q y = n). As one sees 
there are well-pronounced dispersion curves oc q 2 in good 
agreement with experiment [23,24]. Note, that in the 
RPA the dispersion of the resonance is much weaker (in 
particular the upper branch of the dispersion) due to a 
(5-function character of the resonance condition [21] while 
Z— and n-tcrms, in particular, lead to a redistribution of 
the spectral weight away from Q = (tt,it). Note, due to 
the tetragonal symmetry the same dispersion takes place 
for q y (q x = n). 

The position of the resonance peak is determined 
mainly by the magnitude of the rf x 2_ y 2-wave supercon- 
ducting gap, Ao, the superexchange coupling constant, 
J, and by the proximity of the Fermi energy to the ex- 
tended saddle point (the so-called van-Hove singularity) 
in the density of states determined by the ratio of t'/t. 
In order to illustrate this dependence we show in Fig. 
3 the resonance peak position as a function of J/t and 
Ao/i at optimal doping. One can clearly see that it de- 
pends almost linearly on 2 Ao . This follows from our Eq. 
(11). In particular, the superconducting gap determines 
mainly the position of the continuum of the spin excita- 
tions from Imxo an d which is around 2A [20]. Thus, 
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FIG. 2. (a) Calculated frequency and momentum depen- 
dence of Imx(q,o>) away from Q = (n, n). (b) The dispersion 
of the resonance peak. Two branches of the dispersion curves 
are in good agreement with recent experimental data [23,24]. 
Note, the solid curve is a guide to the eye. 

at the fixed value of J/t and t'/t this continuum influ- 
ences also the resonant condition because the difference 
between the position of the resonance peak (SDW col- 
lective mode) and the continuum of the states has to 
remain the same [20]. On the other hand, for the fixed 
values of t'/t and A /i the resonance condition depends 
mainly on the value of J. Furthermore, its dependence 
is more complicated than in the case of Ao. At the rel- 
atively small values of the superexchange coupling the 
resonance lies close to the continuum and is more sensi- 
tive to the change of the superconducting gap magnitude 
than to the superexchange coupling constant. On the 
contrary one could see from the Fig. 3 that if the res- 
onance peak is at small energies and lies relatively far 
from the continuum (which happens for large values of 
J/t) it will be most influenced by the change of J/t val- 
ues. We note that the shift of the resonance peak to the 
lower frequency enhances its intensity and visa versa. In 
addition, the ratio of t'/t which determines the nesting 
of the Fermi surface and the position of van-Hove sin- 
gularity influences also the resonance peak. Its influence 
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FIG. 3. Calculated dependence of the resonance frequency 
uj res on the superexchange coupling constant J (open circles) 
and on the value of the d^^-wave superconducting gap Ao 
(filled squares). Note, the lines are the guides to the eye. 

is somewhat similar to J with one important difference. 
In particular, an increase of t'/t rather weakens the in- 
tensity of the resonance peak than changes the position 
itself (not shown). Note, that our analysis agrees qual- 
itatively with the previous ones [20]. Most importantly 
we find that by the slight variation of all parameters one 
could find another realistic set of A /t, J/t, and t'/t 
which would also correctly explain the position of the 
resonance peak at optimal doping [25]. This means that 
at present stage we could make only a qualitative analysis 
of the experimental data. Further factors like orthorhom- 
bic distortions may influence the difference between the 
resonance peak intensity in I^S^CaCusOs+a; (BSCCO) 
and YBa 2 Cu 3 7 (YBCO). 

Finally, we discuss the influence of the electron-phonon 
interaction on the resonance peak formation by chang- 
ing the isotope mass of 16 O by 18 O. This shifts the 
average frequency of the LO phonon mode and conse- 
quently renormalizes the hopping integral t and the su- 
perexchange coupling constant, J. Most importantly, the 
electron-phonon interaction changes most dramatically 
the hopping integral t rather than the superexchange 
coupling J. In particular, as can be seen from Eq.(6) 
the superexchange coupling constant J changes less than 
1% upon substituting the isotopes [26] which agrees well 
with experimental data [27]. Therefore, there is almost 
no influence of the isotope substitution on the resonance 
peak determined from RPA, since in this approximation 
its formation is determined mainly by J. In particular, 
we find within RPA no change in the oj res value upon 
changing the isotopes. In the case of Eq.(ll) the most 
important contribution to the isotope effect on the reso- 
nance peak appears due to n(q, lu) oc tk- In particular, 
using our estimation given above we find that at opti- 
mal doping the hopping integral changes by 6% upon 
replacing 16 by 18 0. This results in the lowering of the 
resonance frequency at (tt,tt) from 41meV for the 16 
isotope towards 39meV for the 18 O sample. This leads 
to a res = - rf J"^| s ~ 0.4 for optimally-doped cuprates. 
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This effect is beyond the experimental error and can be 
further tested experimentally. Furthermore, in the un- 
derdoped cuprates one may expect larger isotope effect 
due to a larger value of jE* /hui* [13]. At the same time 
the superconducting transition temperature which is de- 
termined by J shows much weaker isotope effect and is 
around «t c ~0.05 [26]. Therefore, even if the supercon- 
ductivity is driven by the magnetic exchange the reso- 
nance peak formation can be significantly renormalizcd 
by the strong electron-phonon interaction. 

To summarize, we analyze the influence of the elec- 
tronic correlations and the electron-phonon interaction 
on the dynamical spin susceptibility in layered cuprates. 
The electronic correlations taken beyond RPA redis- 
tribute the spectral weight of the resonance peak away 
from (-7T, 7r) leading to the pronounced dispersion. This is 
in good agreement with recent INS data [23,24]. Further- 
more, we find the isotope effect on the resonance peak due 
to strong coupling of the carriers to LO phonon mode. 
It results from both electron-phonon coupling and elec- 
tronic correlation effects. In contrast to the small iso- 
tope effect on the superconducting transition tempera- 
ture we find larger isotope coefficient on the resonance 
peak a res « 0.4 in optimally-doped cuprates. We also 
would like to note that the value of the isotope coefficient 
depends strongly on the value of the exponential factor. 
Therefore, the experimental verification of our prediction 
is desirable. In particular, it would put a strong con- 
straint on the ingredients the theory of cuprates must 
contain. 
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